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(57) Abstract 

Methods and systems are described that 
examine tissue (8) positioned between input 
ports (101, 103) and a detection port (107). At 
least one source of a visible or infrared wave- 
length (102, 104) is provided that introduces 
electromagnetic radiation into the subject. The 
detection port is optically coupled to a detec- 
tor (106) that is connected to a detector circuit 
(101). Radiation intensities are selected for in- 
troduction at the input ports to define a null 
plane (105) in the tissue. The detection port is 
positioned relative to the null plane. Radiation 
is introduced into the subject at the first input 
port and the radiation that migrates through the 
tissue is detected. The detector circuit stores 
a first detector signal corresponding to the first 
detected radiation. Radiation is introduced at 
the second input port and is detected. The first 
detector signal is substracted from a second de- 
tector signal corresponding to the second de- 
tected radiation to obtain processed data. 
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OPTICAL TECHNIQUES FOR EXAMIMATION OF BIOLOGICAL TISSUE 

Background of the Invention 
This invention relates to examination and imaging 
5 of biological tissue using visible or infra-red 
radiation • 

Traditionally, potentially harmful ionizing 
radiation (for example , X-ray or y-ray) has been used to 
image biological tissue. This radiation propagates in the 

10 tissue on straight, ballistic tracks, i.e., scattering of 
the radiation is negligible. Thus, imaging is based on 
evaluation of the absorption levels of different tissue 
types. For example, in roentgenography the X-ray film 
contains darker and lighter spots. In more complicated 

15 systems, such as computerized tomography (CT) , a cross- 
sectional picture of human organs is created by 
transmitting X-ray radiation through a section of the 
human body at different angles and by electronically 
detecting the variation in X-ray transmission. The 

20 detected intensity information is digitally stored in a 
computer which reconstructs the X-ray absorption of the 
tissue at a multiplicity of points located in one cross- 
sectional plane. 

Near infra-red radiation (NIR) has been used to 

25 study non-invasively the oxygen metabolism in tissue (for 
example, the brain, finger, or ear lobe) . Using visible, 
NIR and infra-red (IR) radiation for medical imaging 
could bring several advantages. In the NIR or IR range 
the contrast factor between a tumor and a tissue is much 

30 larger than in the X-ray range. In addition, the visible 
to IR radiation is preferred over the X-ray radiation 
since it is non-ionizing; thus, it potentially causes 
fewer side effects. However, with lower energy 
radiation, such as visible or infra-red radiation, the 

35 radiation is strongly scattered and absorbed in 
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biological tissue, and the migration path cannot be 
approximated by a straight line, making inapplicable 
certain aspects of cross-sectional imaging techniques. 

Several different approaches to NIR imaging have 
5 been suggested in the past. One approach undertaken by 
Oda et al. in "Non-Invasive Hemoglobin Oxygenation 
Monitor and Computerized Tomography of NIR Spectrometry," 
SPIE Vol. 1431, p. 284, 1991, utilizes NIR radiation in 
an analogous way to the use of X-ray radiation in an X- 

10 ray CT. In this device, the X-ray source is replaced by 
three laser diodes emitting light in the NIR range. The 
NIR-CT uses a set of photomultipliers to detect the light 
of the three laser diodes transmitted through the imaged 
tissue. The detected data are manipulated by a computer 

15 of the original X-ray CT scanner system in the same way 
as the detected X-ray data would be. 

Different approaches were also suggested by S.R. 
Arriadge et al. in "Reconstruction Methods for Infra-red 
Absorption Imaging," SPIE Vol. 1431, p. 204, 1991; F. A. 

20 Grtinbaum et al. in "Diffuse Tomography," SPIE Vol. 1431, 
p. 232, 1991; B. Chance et al., SPIE Vol. 1431 (1991), p. 
84, p. 180, and p. 264; and others who recognized the 
scattering aspect of the non-ionizing radiation and its 
importance in imaging. None of those techniques have 

25 fully satisfied all needs in tissue examination. 

In summary, there continues to be a need for an 
improved system which utilizes visible or IR radiation of 
wavelengths sensitive to endogenous or exogenous pigments 
to examine or image biological tissue. 

30 Summary of the Invention 

The invention relates to systems and methods for 
spectroscopic examination of a subject positioned between 
input and detection ports of the spectroscopic system 
applied to the subject. 
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According to one aspect, the invention features a 
spectroscopic system for examination of tissue of a 
subject, including; at least one light source of 
electromagnetic radiation of a visible or infrared 
5 wavelength selected to be scattered and absorbed while 
migrating in the tissue; at least two input ports, 
optically coupled to the light source, constructed to 
introduce at selected input locations of the examined 
tissue the radiation of known intensities that define a 

10 null plane in the tissue; a detection port located at a 
selected detection location of the examined tissue 
relative to the null plane; a detector, optically coupled 
to the detection port, constructed to detect during 
operation the radiation that has migrated in the examined 

is tissue; a detector circuit connected to and receiving 
detection signal from the detector; the detector circuit 
including a sample-and-hold circuit and a subtraction 
circuit, both connected to the detector circuit, 
constructed to subtract detection signals corresponding 

20 to radiation that has migrated from a first input port to 
the detection port and from a second input port to the 
detection port, respectively, to obtain processed data; 
and a processor, connected to and receiving the processed 
data from the subtraction circuit, adapted to evaluate 

25 the examined tissue. 

According to another aspect, the invention 
features a spectroscopic system for examination of tissue 
of a subject, including: a source of electromagnetic 
radiation of a visible or infrared wavelength; an input 

30 port, optically coupled to the light source, constructed 
to introduce at a selected input location of the examined 
tissue the radiation; a detector optically coupled to at 
least two detection ports located at selected detection 
locations defining a null plane in the examined tissue, 

35 the detector constructed to detect radiation that has 
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migrated in the examined tissue to the detection ports; a 
detector circuit connected to and receiving detection 
signal from the detector, the detector circuit including 
a sample-and-hold circuit and a subtraction circuit; the 
5 detector circuit constructed to correlate emission of the 
radiation from the input port with detection of radiation 
scattered and absorbed while migrating in the tissue at 
the first detection port, the detected radiation being 
stored as a first detection signal; the detector circuit 
10 further constructed to correlate emission of the 

radiation from the input port with detection of radiation 
scattered and absorbed while migrating in the tissue at 
the second detection port, the detected radiation being 
stored as a second detection signal; the subtraction 
15 circuit constructed to subtract the detection signals; 

and a processor, connected to and receiving the processed 
data from the subtraction circuit, constructed to 
evaluate the examined tissue. 

Embodiments of the invention may include one or 
20 more of the following additional features. 

The spectroscopic system may include intensity 
control means constructed to regulate intensities of 
radiation emitted from the first and second input ports. 
The intensity control means may be constructed to 
25 regulate the intensities in a manner that sweeps the null 
plane over at least a portion of the volume of the 
examined tissue. 

The spectroscopic system may preferably include 
positioning means constructed to displace the detection 
30 port to detection locations corresponding to the null 
plane or positioning means constructed to displace the 
input ports to selected locations. 

The spectroscopic system may preferably include 
detector controller means constructed to changes the 
35 relative sensitivity of detection at the first and second 
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detection port in order to sweep the null plane over at 
least a portion of the volume of the examined tissue. 

Preferably, the subtraction circuit includes an 
analog to digital converter, connected to the sample-and- 
hold circuit, constructed to digitize the detection 
signal to produce digital detection signal, the 
subtraction circuit subtracting the digital detection 
signals corresponding to radiation that has migrated from 
a first input port to the detection port and from a 
second input port to the detection port, respectively, to 
obtain the processed data- The processor may preferably 
be further adapted to locate, in the tissue volume, a 
tissue region exhibiting different scattering or 
absorptive properties than the rest of the examined 

tissue volume. 

The input or detection ports may be preferably 
arranged in a linear array. The input or detection ports 
may be preferably arranged a two dimensional array. The 
spectroscopic may preferably further include an image 
processor, connected to and receiving the processed data 
from the processor, constructed to store processed data 
corresponding to different combinations of input and 
detection ports and create image data; and a display, 
connected to the image processor, constructed to display 
the image data representing the examined tissue. 

The wavelength may be preferably sensitive to an 
endogenous pigment of the examined tissue. The 
wavelength may be preferably sensitive to an exogenous 
pigment of the examined tissue. 

In another general aspect, the invention features 
a scheme for spectroscopic examination of tissue 
including the steps of: providing at least one light 
source of electromagnetic radiation of a visible or 
infrared wavelength selected to be scattered and absorbed 
while migrating in the tissue, the source being optically 
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connected to at least two input ports , and a detection 
port optically connected to a detector, the detector 
connected to a detector circuit; positioning a first 
input port and a second input port relative to selected 
5 input locations of a subject; selecting for each input 
port first and second radiation intensities to be 
introduced to the tissue, the selected radiation 
intensities defining a null plane in the tissue; 
positioning the detection port relative to a selected 
10 detection location of the examined tissue corresponding 
to the null plane, the input locations and detection 
location defining a volume of the examined tissue of the 
subject; introducing into the subject, at the first input 
port, radiation of the first intensity; detecting, at the 
15 detection port, the first radiation that has migrated in 
the examined tissue; storing, in the detector circuit, a 
first detector signal corresponding to the first detected 
radiation; introducing into the subject, at the second 
input port, radiation of the second intensity; detecting, 
20 at the detection port, the second radiation that has 

migrated in the examined tissue; storing, in the detector 
circuit, a second detector signal corresponding to the 
second detected radiation; subtracting the first detector 
signal from the second detector signal to obtain 
25 processed data; and examining the tissue volume using the 
processed data. 

The spectroscopic method may further include the 
step of selecting the first and second radiation 
intensities is preferably performed in a manner that 
30 sweeps the null plane over at least a portion of the 

volume of the examined tissue and the step of positioning 
the detector to detection locations corresponding to the 
swept null plane. 

In another general aspect, the invention features 
35 a scheme for spectroscopic examination of tissue 
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including the steps of: providing a source of 
electromagnetic radiation of a visible or infrared 
wavelength selected to be scattered and absorbed while 
migrating in the tissue, the source being optically 
5 coupled to an input port, and providing at least two 
detection ports optically coupled to at least one 
detector, the detector connected to a detector circuit; 
positioning the input port relative to selected input 
locations of the tissue; positioning a first detection 

10 port and a first detection port relative to selected 

detection locations of the examined tissue, the locations 
defining a null plane in the tissue a volume of the 
examined tissue of the subject; introducing into the 
tissue, at the input port, radiation of a selected 

15 intensity and a selected wavelength; detecting, at the 
first detection port, radiation that has migrated in the 
examined tissue and storing, in the detector circuit, a 
first detector signal corresponding to the detected 
radiation; detecting, at the second detection port, 

20 radiation that has migrated in the examined tissue and 
storing, in the detector circuit, a second detector 
signal corresponding to the detected radiation; and 
subtracting the first detector signal from the second 
detector signal to obtain processed data corresponding to 

25 properties of the tissue volume. 

The detecting steps may be performed in a manner 
that changes the relative sensitivity of detection at the 
first and second detection port in order to sweep the 
null plane over at least a portion of the volume of the 

30 examined tissue and the method may further include, 

simultaneously with the sweeping, positioning the input 
port to input locations corresponding to the swept null 
plane. 

Further embodiments of the invention may include 
35 one or more of the following features. 
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The first detector signal and the second detector 
signal are preferably stored and subtracted in an analog 
form by the detection circuit • The spectroscopic method 
preferably further includes, before the storing steps, 
5 converting the first and second detection signals to a 
digital form, the subtracting step being performed 
digitally on the first detector signal and the second 
detector signal in an digital circuit. 

The steps of positioning the input ports and the 

10 detection port at the respective selected locations 
preferably include placing the ports on the surface of 
the examined tissue. The steps of positioning the input 
ports at the respective selected locations preferably 
include orienting the input ports relative to the input 

15 locations thereby enabling introduction of the radiation 
at the input locations of the examined tissue. The step 
of orienting the input ports preferably include directing 
an optical element of the input ports to the input 
locations. The positioning steps are preferably 

20 performed in a manner that sweeps the null plane and the 
detector over at least a portion of the volume of the 
examined tissue. 

The spectroscopic method may further include 
locating, in the tissue volume, a tissue region 

25 exhibiting different scattering or absorptive properties 
than the rest of the examined tissue volume. The method 
may preferably further include imaging the examined 
tissue including the tissue region of different 
scattering or absorptive properties. The method may 

30 preferably further include displaying an image of the 
examined tissue by utilizing to the processed data and 
relative locations of the input ports and the output 
ports . 

The method may preferably further include a step 
35 of introducing an exogenous pigment into the tissue and 
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selecting a wavelength being sensitive to the pigment. 
The exogenous pigment may be preferentially accumulated 
in a tissue region exhibiting different scattering or 
absorptive properties. The exogenous pigment may be 
5 fluorescing when irradiated by selected wavelength, and 
the detecting steps may detect preferentially radiation 
wavelength of the fluorescing pigment. 

In general, according to another aspect of the 
invention, a spectroscopic system includes at least one 

10 light source adapted to introduce, at multiple input 

ports, electromagnetic non-ionizing radiation of a known 
time-varying pattern of photon density and of a 
wavelength selected to be scattered and absorbed while 
migrating in the subject, the input ports being placed at 

15 selected locations on the subject to probe a selected 
quality of the subject; radiation pattern controller 
adapted to achieve selected a time relationship of the 
introduced patterns to form resulting radiation that 
possesses a substantial gradient in photon density as a 

20 result of the interaction of the introduced patterns 
emanating from the input ports, the radiation being 
scattered and absorbed in migration paths in the subject. 
The system also includes a detector adapted to detect 
over time, at a detection port placed at a selected 

25 location on the subject, the radiation that has migrated 
in the subject; processor adapted to process signals of 
the detected radiation in relation to the introduced 
radiation to create processed data indicative of the 
influence of the subject upon the gradient of photon 

30 density; and the processor (evaluation means) adapted to 
examine the subject by correlating the processed data 
with the locations of the input and output ports. 

According to another aspect of the invention, a 
spectroscopic system includes at least one light source 

35 adapted to introduce, at multiple input ports, 
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electromagnetic non-ionizing radiation of a known time- 
varying pattern of photon density and of a wavelength 
selected to be scattered and absorbed while migrating in 
the subject, the input ports being placed at selected 
5 locations on the subject to probe a selected quality of 
the subject; radiation pattern controller adapted to 
achieve a selected time relationship of the introduced 
patterns to form resulting radiation that possesses a 
substantial gradient in photon density as a result of the 

10 interaction of the introduced patterns emanating from the 
input ports, the radiation being scattered and absorbed 
in migration paths in the subject. The system also 
includes a detector adapted to detect over time, at a 
detection port placed at a selected location on the 

15 subject, the radiation that has migrated in the subject; 
displacement means adapted to move the detection port to 
various locations on a predetermined geometric pattern, 
the various locations being used to detect over time 
radiation that has migrated in the subject; processor 

20 adapted to process signals of the detected radiation in 
relation to the introduced radiation to create processed 
data indicative of the influence of the subject upon the 
gradient of photon density; and the processor (evaluation 
means) adapted to examine the subject by correlating the 

25 processed data with the locations of the input and output 
ports . 

According to another aspect of the invention, a 
spectroscopic system includes at least one light source 
adapted to introduce, at multiple input ports, 

30 electromagnetic non-ionizing radiation of a known time- 
varying pattern of photon density and of a wavelength 
selected to be scattered and absorbed while migrating in 
the subject, the input ports being placed at selected 
locations on the subject to probe a selected quality of 

35 the subject; radiation pattern controller adapted to 
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achieve a selected time relationship of the introduced 
patterns to form resulting radiation that possesses a 
substantial gradient in photon density as a result of the 
interaction of the introduced patterns emanating from the 
5 input ports, the radiation being scattered and absorbed 
in migration paths in the subject. The system also 
includes at least one detector adapted to detect over 
time, at multiple detection ports placed at selected 
locations on the subject, the radiation that has migrated 

10 in the subject; processor adapted to process signals of 
the detected radiation in relation to the introduced 
radiation to create processed data indicative of the 
influence of the subject upon the gradient of photon 
density, and the processor (evaluation means) adapted to 

15 examine the subject by correlating the processed data 
with the locations of the input and output ports. 

Preferred embodiments of this aspect of the 
invention include displacement means adapted to move at 
least one of the detection ports to another location on a 

20 predetermined geometric pattern, the other location being 
used to perform the examination of the subject. 

Preferred embodiments of this aspect of the 
invention include rotation means adapted to rotate 
synchronously the optical input ports while introducing 

25 the resulting radiation along a predetermined geometric 
pattern, the input port rotation being used to perform 
the examination of a region of the subject. 

According to another aspect of the invention, a 
spectroscopic system includes a light source adapted to 

30 introduce, at an input port, electromagnetic non-ionizing 
radiation of a known time-varying pattern of photon 
density and of a wavelength selected to be scattered and 
absorbed while migrating in the subject, the input port 
being placed at a selected location on the subject to 

35 probe a selected quality of the subject; detectors 
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adapted to detect over time, at multiple detection ports 
placed at selected locations on the subject, the 
radiation that has migrated in the subject; the time 
relationship of the detection over time, at the detection 
5 ports, being selected to observe a gradient in photon 
density formed as a result of the interaction of the 
introduced radiation with the subject. The system also 
includes processor adapted to process signals of the 
detected radiation in relation to the introduced 

10 radiation to create processed data indicative of the 
influence of the subject upon the gradient of photon 
density, and the processor (evaluation means) adapted to 
examine the subject by correlating the processed data 
with the locations of the input and output ports. 

is According to another aspect of the invention, a 

spectroscopic system includes a light source adapted to 
introduce, at an input port, electromagnetic non-ionizing 
radiation of a known time-varying pattern of photon 
density and of a wavelength selected to be scattered and 

20 absorbed by a fluorescent constituent while migrating in 
the subject, the input port being placed at a selected 
location on the subject to locate the fluorescent 
constituent of the subject; detectors adapted to detect 
over time, at multiple detection ports placed at selected 

25 locations on the subject, fluorescent radiation that has 
migrated in the subject. The system also includes 
processor adapted to process signals of the detected 
radiation in relation to the introduced radiation to 
create processed data indicative of location of the 

30 fluorescent constituent of the subject, and the processor 
(evaluation means) adapted to examine the subject by 
correlating the processed data with the locations of the 
input and output ports. 

In certain preferred embodiments, the 

35 spectroscopic system further includes an image processor, 
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connected to and receiving the processed data from the 
processor, constructed to store processed data 
corresponding to different combinations of input and 
detection ports and create image data, the image data 
5 including data of the tissue region; and a display, 

connected to the image processor, constructed to display 
the image data representing the examined tissue including 



A displacement mechanism is adapted to move 
10 synchronously the optical ports and the detection ports 
to another location on a predetermined geometric pattern; 
this other location is used to perform the examination of 
the subject. 

The spectroscopic system also uses a wavelength 
15 sensitive to endogenous or exogenous pigments of the 
examined biological tissue. 

The spectroscopic system also used to locate a 
fluorescent constituent of interest in the subject; the 
wavelength of the introduced radiation is selected to be 
20 absorbed in the fluorescent constituent, the detected 

radiation is emitted from the fluorescent constituent and 
processed to determine location of the fluorescent 
constituent . 

The time-varying pattern of resulting radiation is 
25 formed by the intensity modulated radiation introduced 
from each of the input ports having selected phase 
relationship that produces in at least one direction a 
steep phase change and a sharp minimum in the intensity 

of the radiation. 
30 The phase relationship of the introduced radiation 

patterns is 180 degrees. 

The modulation frequency of the introduced 
radiation has a value that enables resolution of the 
phase shift that originates during migration of photons 
35 in the subject. 
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Other features and advantages will become apparent 
from the following description and from the claims* 



Brief Description of the Drawing 
Figs. 1 and 1A show diagrammatically a phase 
5 modulation imaging system including several input ports 
and one detection port in accordance with the present 
invention. 

Fig. 2 is a block diagram of the phase modulation 
imaging system including several input ports and several 
10 detection ports in accordance with the present invention. 

Fig. 2A depicts a phased array antenna that 
radiates a directional beam. 

Fig. 2B depicts sequencing of the phases of an 
antiphase multi-element array to achieve an electronic 
15 scan of the photon density gradient in accordance with 
the present invention. 

Fig. 2C depicts four element antiphased array 
designed for a conical scan of the photon density 
gradient in accordance with the present invention. 
20 Fig. 2D depicts the input and output port 

arrangement of an imaging system in accordance with the 
present invention. 

Fig. 3 depicts a phase modulation imaging system 
including an input port and several array detection ports 
25 in accordance with the present invention. 

Fig. 4 is a block diagram of an alternative 
embodiment of a dual wavelength PMS system. 

Fig. 4 A is a schematic diagram of an oscillator 
circuit of Fig. 4. 
30 Fig. 4B is a schematic diagram of a PMT heterodyne 

modulation and mixing network shown in Fig. 4. 

Fig. 4C is a schematic diagram of an AGC circuit 
shown in Fig. 4. 
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Fig. 4D is a schematic diagram of a phase detector 

circuit shown in Fig. 4. 

Fig. 5A, 5B, and 5C illustrate changes in optical 
field propagating in a strongly scattering medium which 
5 includes a strongly absorbing component. 

Fig. 6 shows an experimental arrangement of a two 
element phased array used in an interference experiment. 

Figs. 6A, 6B, and 6C show detected interference 
patterns of two diffusive waves. 
10 Fig. 7 displays the phase shifts measured for a 

two element array (curve A) , and for a single source 
(curve B) . 

Fig. 8A depicts an experimental arrangement of 
sources of a four element phased array and a detector. 
15 Figs. 8B and 8C display the intensities and the 

phase shifts measured for the four element array of Fig. 

8A, respectively. 

Fig. 9A depicts an experimental arrangement of 
sources of a four element phased array, a detector, and a 

20 strongly absorbing object. 

Fig. 9B, 9C display respectively the intensities 
and the phase shifts measured for the four element array 
of Fig. 9A scanning absorbing objects of different sizes. 

Fig. 9D displays the phase shifts measured for the 
25 four element array of Fig. 9 A scanning absorbing objects 
of different absorption coefficients. 

Fig. 10A displays an experimental arrangement of 
sources of a four element phased array, a detector, and 
two strongly absorbing objects. 
30 Fig. 10B displays the phase shifts measured for 

the four element array of Fig. 10A scanning two absorbing 
objects of different sizes. 

Fig. 11 shows diagrammatically a low frequency 
imaging system utilizing a one dimensional source array 
35 in accordance with the present invention. 
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Fig* 11A shows a circuit diagram of the low 
frequency imaging system of Fig* 11. 

Figs. 12 and 13 show diagrammatically a low 
frequency imaging systems utilizing a two dimensional 
5 source array in accordance with the present invention. 

Figs. 14 and 14A show a scanning system 
constructed for imaging of breast tissue. 



Description of the Preferred Embodiments 
Imaging system embodiments of the present 

10 invention based upon interference effects of radiation 
migrating in a subject having scattering and absorptive 
properties are shown in Figs. 1, 2, and 3. The systems 
effectively utilize, in this scattering medium, a 
directional beam of visible or IR radiation generated 

15 and/or detected by an array of sources and/or detectors, 
respectively. For instance, in the case of an array of 
sources, each source is placed at a selected location in 
the array and emits intensity modulated radiation, 
preferably coherent radiation from a laser diode, of a 

20 selected intensity and phase. The criteria for selecting 
the source locations, the intensities, and the phases of 
the respective sources is the shape of the desired beam 
that at any time point possesses a substantial photon 
density gradient produced by interference effects of 

25 radiation from the various sources. This gradient of 
photon density is localized and has directional 
properties. Overall, the resulting radiation formed by 
interference of the radiation of the individual sources 
migrates in a selected direction in the subject. In an 

30 antiphase system, the wavefront of the beam has sections 
of equal photon density separated by a sharp localized 
change in photon density. Selected different locations 
of the photon density gradient are shown in Fig. 2B. 



WO 96/16596 



PCT/US95/15694 



- 17 - 

In general, the wavefront propagates in the 
selected direction in the subject and the gradient of 
photon density is localized in one or more planes 
extending from the source array in a selected direction, 
5 If the subject includes a localized object having 

different scattering and absorptive properties from those 
of the surrounding environment, the propagating radiated 
field is perturbed. This perturbation is detected and 
from the source detector geometry the perturbing object 

10 can be located. 

In one preferred embodiment shown in Figs. 1 and 
1A, the imaging system utilizes an array of laser diodes 
12, 14, 16, and 18 for introducing light into the tissue 
at selected locations. The geometry of optical input 

15 ports 11, 13, 15, 17 and of an optical output port 19 is 
selected to examine a specific part of the tissue. From 
the known geometry of the optical input ports and the 
detection port and from the shape of the introduced and 
detected radiation, a computer can locate a hidden object 

20 9 of examined tissue 8 (For example, the head or breast) . 
A master oscillator 22, which operates at 200 MHz, 
excites laser diode 12 through 18, that emit light of a 
selected wavelength (e.g., 760 nm) . The light from each 
laser diode is conducted to the respective input port 

25 placed on a subject via a set optical fibers. A detector 
24 detects the light that has migrated through the 
examined tissue. Preferably, detector 24 includes a 
photomultiplier tube (e.g., Hamamatsu R928) powered by a 
high voltage supply which outputs about 900 V in order to 

30 ensure a high gain. A local oscillator 26 operating at a 
convenient offset frequency (e.g., 25 KHz) sends a signal 
to a mixer 28 and a reference signal to detector 24. 
Accordingly, an output waveform 25 from detector 24 is at 
a carrier frequency equal to the difference of the 

35 detected and reference frequency, i.e., 25 KHz. 
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Detector 24 (for example, PMT Hamamatsu R928 or 
Hamamatsu R1645u) detects the scattered and absorbed 
light that has migrated through the subject. Detection 
port 19 is located several centimeters from the location 
5 of the input ports. The PMT detector is connected to the 
subject by the fiber optic guide, or, alternatively, may 
be directly placed on the subject. It has been found 
that the most cost-effective detector for measuring 
signals of frequencies on the order of 10 8 Hz is 

10 Hamamatsu R928. However, the Hamamatsu R1645u detector 
is preferred due to its high precision. The second 
dynode of the PMT of detector 24 is modulated by 200.025 
MHz signal 27 so that the 25 KHz hetrodyned signal 25 is 
received by a phase detector 30. Phase detector 30 also 

15 receives reference signal 29 from mixer 28. If phase 
detector 30 is a lock-in amplifier then the output 
signals are the phase shift and the intensity of the 
detected signal. Both the phase shift and the intensity 
of the detected light characterize the migration path of 

20 photons in the subject (e.g., the brain tissue). 

Alternatively, a tunable dye laser or other laser 
source connected to a wide band acousto-optical modulator 
operating at the carrier frequency, e.g., 200 MHz can be 
used instead of the laser diode. The acousto-optical 

25 modulator modulates the intensity of the light emitted by 
the laser at the selected carrier frequency. 

The invention also envisions using only one source 
of coherent light that irradiates one end of several 
optical fibers at the same time. The other end of each 

30 fiber is placed on the subject at a selected input port 
location. This source radiates light of a selected time 
varying pattern. The phase relationship and the 
intensity of the light carried by each fiber is varied by 
creating a time delay (e.g., different fiber length) and 

35 by coupling different amounts of light into each fiber. 
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The imaging systems of Figs. 1, 2, and 3 are shown 
to have a light source of a single wavelength; however, a 
dual wavelength imaging system is also envisioned 
according to this invention. In the dual wavelength 
5 imaging system two laser diodes or a tunable wavelength 
laser generate light of two wavelengths that is coupled 
to an optical fiber. Such a system will now be 
described. 

A dual wavelength operation is shown in Fig 4. 

10 The system includes a master oscillator 60 operating at 
200 MHz and an oscillator 62 operating at 200.025 MHz 
which is offset 25 KHz from the master oscillator 
frequency. The offset frequency of 25 KHz is a 
convenient frequency for phase detection in this system; 

15 however, other offset frequencies as high as a few 
megahertz can be used. Oscillator 60 alternatively 
drives two laser diodes 64 and 66 using switches 61, 
61a,... These switches are driven electronically to 
couple a selected wavelength into the optical fiber and 

20 also to achieve a selected radiation pattern resulting 
from the radiation emanating from the individual fibers. 
An output 8 mm fiber coupler 72 collects photons for an 
R9 2 8 PMT detector 7 4 . The second dynode ( shown in Fig . 
3B) of PMT 74 is modulated with a 200.025 MHz reference 

25 signal generated by oscillator 62 and amplified by an 
amplifier 63. Thus, the output signal of the PMT 
detector has a frequency of 25 KHz. PMT detector 74 
alternately detects light of the two laser diodes that 
has migrated in the tissue and produces corresponding 

30 output signals, which are filtered by a filter 78 and 

leveled by an automatic gain control (AGC) circuit 79. A 
reference signal of 25 KHz is produced in a mixer 65 by 
mixing the 200 and 200.025 MHz oscillator signals. The 
reference 25 kHz signal is also leveled using the second 

35 AGC 77 and fed into a phase detector 80. Phase detector 
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80 generates a signal indicative of the phase of each 
output signal relative to the phase of the reference 
signal. The outputs of phase detector 80 are alternately 
selected by an electronic switch 82, filtered, and then 
5 input to an adder 84 and a subtractor 86 to produce sum 
and difference signals proportional to <L>xi +<L> X2 and 
<L>Xi~ <Ij> X2* T ^e difference and sum signals are then used 
to calculate changes in the probed pigment and in the 
blood volume, respectively. 

io A schematic diagram of preferred oscillator 60 or 

62 is shown in Fig 4A. This circuit has a drift of only 
0.03 degrees/hr. (Weng, et al., "Measurement of 
Biological Tissue Metabolism Using Phase Modulation 
Spectroscopic Measurement," SPIE, Vol. 143, p. 161, 1991, 

15 which is incorporated herein by reference) . The crystal 
is neutralized, which enables operation at resonance, and 
thus achieves long-term stability. The respective 
crystals of oscillators 60 and 62 are offset from each 
other by 25 kHz. This circuit provides a sufficient 

20 output to directly drive a 5 mW laser diode. 

A modulation circuit 75 for the second dynode of 
the PMT is shown in Fig. 4B. This circuit uses a 
resonant circuit 75a with an impedance of 20,000 ohms 
instead of the usual 50 n load with very high power 

25 dissipation, providing a 50 V drive of the 

photomultiplier dynode while dissipating only a few watts 
of power. 

To obtain stable operation of the phase detector, 
a stable input signal is required. The 25 KHz AGC 

30 circuit 77 illustrated in Fig. 4C includes an MC 1350 
integrated circuit Ul, featuring wide range AGC for use 
as an amplifier. The signal amplitude is controlled by a 
feedback network, as shown. A major reason for the 
accurate detection of phase changes by the PMT system is 

35 that the phase detector input signal level is kept nearly 
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constant by the AGC circuit. Since the input voltage 
change of between 2 and 6 volts causes variation in the 
phase shift of only 0.2%, the AGC circuit eliminates the 
need for a very stable high voltage power supply. 
5 A preferred phase detector circuit is shown in 

Fig. 4D. Two sinusoidal signals (the measurement signal 
and the reference signal) are transformed to a square 
wave signal by a Schmitt trigger circuit 80a. The phase 
of the square wave signal is shifted by an RC change 

10 (composed of Rll, R12, C8) , which makes it possible to 
change the measuring range. The detector further 
includes a 74HC221 integrated circuit. The lock- in 
amplifier technique obtained to derive the difference of 
the phase and amplitude of the two signals has the 

15 highest signal to noise ratio possible for this type of 
equipment. 

The above-described systems utilize the carrier 
frequency on the order of 10 8 Hz which is sufficiently 
fast to resolve the phase shift of the detected light. 

20 The characteristic time, the time it takes for a photon 
to migrate between an input port and an output port, is 
several nanoseconds. The sensitivity of the system is 
high, approximately 70° per nanosecond or 3° per 
centimeter change of pathlength, as observed in 

25 experimental models. Selection of the modulation 

frequency also depends on the desired penetration depth 
and resolution of the imaging system that will be 
described below. If deep penetration is desired, low 
modulation frequency (e.g., 40 MHz) is selected, and if 

30 shallow penetration is needed, modulation frequencies as 
high as 10 9 Hz can be used. 

Referring to Figs. 1 and 1A, a master oscillator 
22 operates at a modulation frequency in the range of 40 
to 4 00 MHz selected according to the desired penetration 

35 depth of the optical field. The array of laser diodes 
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12 , 14, 16, and 18 generates a highly directional 
radiation pattern, which is employed in the tissue 
examination . 

In one preferred mode of operation, laser diodes 
5 12 to 18 operate in a phased array pattern which is 

introduced into the tissue and detected by a single PMT 
detector 30. Master oscillator 22 operating at 200 MHz 
drives a multi-channel phased splitter which gives 
outputs at predetermined phases. Input ports 11 through 

10 17 are located at selected distances and an appropriate 
phasing of the array creates a directional beam and 
enables scanning of the optical field in two dimensions 
across the tissue, as shown in Figs. 2A, 2B, and 2D. 
After migrating through the tissue, the optical field is 

15 collected in a large area fiber on selected locations 19. 
The detected signals are heterodyned in the PMT detector 
24 by utilizing the output of local oscillator 26, 
operating at a 25 kHz offset frequency, to detector 24. 
The resulting 25 kHz signal is phase detected with 

20 respect to the output signal 29 of mixer 28 and detector 
24. Phase detector 30 outputs the phase and the 
intensity of signal 25. The detected phase shifts and 
intensities are stored and used for construction of an 
image of the subject. This is performed by computer 

25 control 34, which governs the operation of the system. 

Fig. 2 depicts a phase modulation imaging system 
comprising an input port array for introducing radiation 
and detection port array for detecting radiation that has 
migrated in the subject. The operation of the system is 

30 controlled by computer control 34, which coordinates a 
transmitter unit 32 with a receiver unit 42. Transmitter 
unit 32 comprises several sources of visible or IR 
radiation adapted to introduce a selected time-varying 
pattern of photon density into subject 8 by array of 

35 input ports 31, 33, 35, and 37. Receiver unit 4 2 detects 
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radiation that has migrated in the subject from the input 
port array to an array of detectors 39 , 41, 42, and 47. 

The radiation sources of transmitter unit 32 are 
intensity modulated at a frequency in the range of 40 MHz 
5 to 200 MHz, as described for the imaging system of Fig. 
1. Receiver unit 42 detects and processes the radiation 
using the same principles of the phase and amplitude 
detection as described above. The signal detected at 
individual ports can be phased using appropriate delays. 

io Several modes of operation of the transmitter 

array and receiver array are described in Figs. 2A, 2B, 
2C, and 2D. Referring to Fig. 2A, it has been known, 
that for a simple horizontal linear array of N identical 
elements radiating amplitude modulated light spaced a 

15 distance, d, apart. The radiating wavefront is created 
by the interference effect. If all elements radiate in 
phase the wavefront propagates in a direction 
perpendicular to the array. However, by appropriately 
phasing the radiating elements, the resulting beam can 

20 scan space in two dimensions. We consider the phases of 
the signal along the plane A-A whose normal makes an 
angle 6 0 with respect to the array normal. The phase of 
the signal from the first radiator lags the phase of the 
second radiator by a phase angle (27r/X)d sin 8 0 because 

25 the signal from the second radiator has to travel a 



d sin 6 0 longer than the signal from the first radiator 
to reach plane A-A. Similarly, the phase of the signal 
from the n th radiator leads that from the first radiator 

30 by an angle n(27r/X))d sin 6 0 . Thus, the signals from the 
various radiators can be adjusted to be in-phase along 
the A-A plane, if the phase of each radiator is increased 
by (27r/A)d sin 6 0 . Consequently, at a point on the 
wavefront in the far field of the antenna array the 

35 signals from the N radiators will add up in phase, i.e., 
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the intensity of the total normalized signal is a sum of 
the signals from the individual sources. The constructed 
pattern has a well defined directional characteristic and 
a well pronounced angular dependence, i.e., the antenna 
5 pattern has a well defined transfer characteristic of the 
antenna with respect to the angle 6 0 . 

Fig. 2B depicts an arrangement of phases for the 
sources the system of Fig. 2 operating in one preferred 
mode of operation. The array of five sources is divided 

10 into two or more portions that are phased 180° apart. 
Each portion has at least one source. The sources of 
each portion radiate amplitude modulated light of equal 
intensity and are spaced so that the resulting beam of 
two or more equally phased sources has a substantially 

15 flat wavefront, i.e., no gradient of photon density. On 
the other hand, there is a sharp 180° phase transition, a 
large gradient in photon density between two antiphased 
portions of the array. Thus, the radiated field 
possesses an amplitude null and a phase transition of 

20 180°, which is due to the large gradient of photon 
density. 

Electronic scanning is performed by appropriately 
varying the apportionment of 0° and 180° phases on the 
sources. The five element array of Fig. 2B can have the 

25 180° phase transition along four different parallel 

planes extending from the array. Scanning is achieved by 
electronically switching the sources by 180°, so that the 
photon density gradient moves in the direction parallel 
to the location of the sources. 

30 Using the principles described in Figs. 2A and 2B, 

a conical scan of a directional beam possessing at least 
one substantial photon density gradient can be 
accomplished using a four element antiphased array, as 
shown in Fig. 2C. The laser diodes are antiphased using a 

35 push pull transformer. The phasing and amplitude of four 



WO 96/16596 



PCT/US95/15694 



laser diodes S lf S 2 , S 3 , and S 4 arranged into a two 
dimensional array is modified sequentially using the 
switches Sw x , Sw 2 , Sw 3 , and Sw 6 and inductances L 1# L 2 , 
L3 , and • 

5 Fig. 2D shows a possible arrangement of the 

transmitter array and the receiver array. The above 
described directional beam enters subject 8 at the 
transmitter array location and is pointed to hidden 
absorber 9 which perturbs the migrating beam. The field 

10 perturbation is measured by the receiver array. Scanning 
of the transmitter array or the receiver array is 
envisioned by the present invention. 

A hidden absorber that includes a fluorescent 
constituent is detected using a selected excitation 

15 wavelength of the laser sources of the transmitter array. 
Then, the radiation is absorbed, and almost instantly a 
fluorescent radiation of a different wavelength is re- 
emitted. The re-emitted radiation propagating in all 
directions is detected by the receiver array. 

20 Fig. 3 depicts a phase modulation imaging system 

comprising one input port and several arrays of detection 
ports. This system operates comparably to the systems of 
Figs. 1 and 2. The 754nm light of a laser diode 48 is 
amplitude modulated using master oscillator 22. The 

25 light is coupled to subject 8 using an input port 49. 

The amplitude modulated light migrates in the subject and 
is scattered from hidden object 9. It is also expected 
that hidden object 9 has a different effective index of 
refraction than subject 8. The migrating radiation is 

30 governed by the laws of diffusional wave optics that are 
described below. The scattered radiation migrates in 
several directions and is detected by detection systems 

50 , 52, and 54. 

Ports 51, 53, and 55 of the detection systems 
35 include either large area fibers or arrays if detection 
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ports. If large area fibers are used then detector 
systems 50, 52, and 54 correspond to detector 24 of Fig. 
1. If arrays detection ports are used, then each of 
detector systems 50, 52, and 54 includes several 
5 individual PMT detectors. The PMT detectors of each 
detector system are phased utilizing a selected phase 
mode, as described above. The phasing is controlled by 
the computer control. The detected signals are 
heterodyned at the PMT's and sent to a phase detector 58. 
10 Phase detector 58 detects alternatively the heterodyned 
signals using a switch 56. Operation of phase detector 
58 is similar to the operation of phase detector 3 0 of 
Fig. 1. The detected phase and amplitude are 
alternatively sent to the computer control using a switch 
15 56a. Even thought only one phase detector is shown in 
Fig 3, the invention envisions use of several phase 
detectors . 

If hidden absorber 9 includes a fluorescent 
constituent, laser diode 4 8 is selected to introduce an 
excitation wavelength (e.g., 754nm) . The introduced, 
intensity modulated radiation, excites the fluorescent 
constituent which re-emits radiation in all directions, 
as shown in Fig. 3. The re-emitted radiation is detected 
using detector systems 50, 52, and 54. To increase the 
system resolution, each detector can be furnished with an 
interference filter selected to pass only the fluorescent 
radiation. 

The interference of several waves, as described in 
Fig 2A, has been long known in a non-scattering medium, 
wherein the radiation propagates on a straight line, but 
not in a strongly scattering medium* Referring to Figs. 
6, 6A, 6B, and 6C, in a simple experiment, interference 
of two different diffusive waves in a strongly scattering 
medium was demonstrated. Propagation of visible IR 
radiation in a scattering medium such as tissue can be 
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described by diffusion of photons, and thus we describe 
it as a diffusive wave. 

Referring to Fig. 6, the two laser diodes were 
separated at a distance of 4 cm and 1.2 cm from the 
s detection port. The intensity modulated light of the two 
laser diodes at freguency 200 MHz was sent through two 
optical fibers to a container with an Intralipid" 
suspension. The source detector distance was varied by 
moving the optical port of the detection fiber along a 

10 line parallel to the position of the sources. Figs. 6A, 
6B, and 6C show measured maxima and minima of the optical 
field migrating in the medium. This data demonstrates 
interference between two diffusive waves created by two 
coherent emitting sources of phase difference 180 

is degrees. Fig. 7 summarizes the experiment, wherein the 
displacement of the detector is plotted against the phase 
shift measured by the detector. The phase shift displays 
the steepest part of the trace, curve A, (slope of about 
360°/cm) at the displacement of about 2.25 cm. Curve B 

20 is measured with an optical field of source S 2 . Here, the 
measured slope is about 30° /cm. When comparing curves A 
and B we demonstrate much higher sensitivity of the null 
detection of the two element array contrasted with a 
diminished sensitivity to the detector displacement when 

25 using a single source arrangement. The sensitivity of 
the two source arrangement is increased by about a factor 
of 10. The sensitivity is further increased when using 
four or more element phased array, which sharpens the 
photon density gradient and thus provides a higher 

30 resolution for locating the hidden object. 

In a strongly scattering medium, the emitted 
photons undergo a large number of collisions and their 
migration can be determined by applying the diffusion 
eguation. The diffusion eguation for photons in a 

35 uniformly scattering medium was solved by E. Gratton et 
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al., "The possibility of a near infrared optical imaging 
system using frequency domain methods." in Mind Brian 
Imaging Program, Japan 1990; and by J. Fishkin et al., 
"Diffusion of intensity modulated near-infrared light in 
5 turbid media", SPIE Vol. 1413 (1991) p. 122. A solution 
of the diffusion equation was obtained for the light of a 
point source (at r=0) radiating S ( 1 + M exp[- 

i(o>t+«)] photons, wherein S is the source strength 
( photons/ sec. ) , M is the modulation of the source at 
xo frequency g>, and e is an arbitrary phase. The photon 
intensity can be calculated as 

I(£,t)= c * p(r,t) , 

wherein p(r,t) is the photon density and c = 10 8 m/s is 
the velocity of light. 
15 When solving the diffusion equation using a 

spherical-harmonics approximation in a non-absorbing 
medium for the density of photons p(r,t) than 

I(r,t) = (I 0 /Dr) + (I 0 /Dr)exp[-r(ca/2cD) , «] 

x exp[ir ((i)/2cD) 3j - i(ut+«)]# 

20 wherein the diffusion constant D is 14 of the mean free 
path. In the absence of an amplitude modulated signal 
(<d=0) the solution corresponds to a spherical wave 
propagating without attenuation. For a non-zero 
frequency, the amplitude of the signal at a frequency <*) 

25 decreases exponentially. The light wave front the 
emitted advances at the constant velocity V 

V = (2Dco) h 



and has wavelength 

30 

X = 2jt(2cD/g>) 3 * 
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The above equations show that higher modulation 
frequencies yield shorter effective wavelengths, and 
smaller diffusion constants also give shorter effective 
wavelengths. In principle, short wavelengths can be 
5 obtained using high frequency modulated waves in a very 
turbid medium. However, the amplitude of the modulated 
wave decreases exponentially with the modulation 
frequency. Therefore, the best resolution, i.e., the 
shortest wavelength, is obtained using the highest 

10 frequency which still gives a measurable signal. The 
diffusion process limits the penetration depth at any 
given modulation frequency, because of the exponential 
decrease of the wave's amplitude, and also decreases the 
velocity of light propagation. 

is The above described diffusion wave approach treats 

amplitude modulated light waves in scattering media using 
the framework of wave optics. The photon intensity, 
calculated as superposition of different waves, 
constitutes a scalar field, propagating at a constant 

20 velocity. At any given modulation frequency, the wave 
optics phenomenology of scalar fields is valid. 
Therefore, in the frequency-domain, the measurement and 
analysis of light diffusing in tissues from several 
sources will undergo constructive and destructive 

25 interference. Furthermore, wave refraction occurs at a 
boundary between two different tissues. It causes a 
deviation of the direction of propagation of the wave 
front, and thus there is a change in the amplitude and 
phase shift of the propagation wave. The direction 

30 change is a function of the ratio of the effective index 
of refraction in the two tissues. In diffusional wave 
optics, on the other hand, the wave's amplitude is 
exponentially attenuated as the wave propagates in the 
scattering medium. This attenuation is in addition to 
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the exponential attenuation caused by finite absorption 
of the medium. 

Amplitude modulated waves propagate coherently in 
the scattering medium; this is crucial for image 
5 reconstruction. It is possible to accurately measure in 
real time, the average intensity, amplitude, and phase of 
the wave front over a large area using a single detector 
or an array of detectors applying well-established 
frequency-domain methods. 

io The emitters are varied sequentially in phase 

starting with the first emitter in the line and followed 
by subsequent emitters. Each emitter emits a spherical 
wave and propagation of the resultant beam is 
perpendicular to the wavefront. If all the transmitter 

15 delays are equal, the beam travels straight ahead. Delay 
lines which produce variable transmitter delays can be 
used to obtain appropriate phasing for steering the beam 
across the tissue. The same principle can apply during 
reception. 

20 There are two important aspects of imaging as 

envisioned by the present invention. The first is a 
geometrical aspect and the second is phasing of the 
transmitters and receivers. 

It is also possible to construct a two-dimensional 

25 array for two-dimensional pointing (e.g., Fig. 2C) . The 
multiplexing switches used with these arrays can be 
constructed as an integral part of the array and can 
consist of field effect transistors arranged so that 
access to any element may be obtained by the application 

30 of two adverse signals. 

In addition to electronic scanning, the two- 
dimensional scanning can be achieved by moving the array 
of sources and detectors in a regular pre-determined 
pattern in a plane parallel to that being investigated in 

35 the subject. For maximum detection, the detector is 
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places in the plane of the photon density gradient of the 
resulting field created by the array of sources. The 
plane of the photon density gradient is swept as the 
array moves. In this sweeping action, as a strongly or 
5 weakly absorbing object enters the radiation field, the 
detector registers a field imbalance due to the above 
described refraction of the propagating radiation. A 
two-dimensional image is formed by storing the 
information while the probe is moved across the subject. 

10 Several scans in different imaging planes are envisioned 
by the invention. If the system is duplicated or time 
shared in two other faces of a cube, an algorithm would 
be used to provide a 3 -dimensional picture of the object 
by triangulation, as known in the art. The data storage 

15 is accomplished electronically. 

The detector detects the intensity and the phase 
shift of the radiation that has migrated in the subject. 
The phase shift depends on the tissue properties, i.e., 
absorption and scattering. For the low frequencies the 

20 phase shift is proportional to ( (1-g) M s /M a ) ** and for the 
high frequencies proportional to l//x a * To obtain desired 
penetration depth, appropriate frequency for both master 
oscillator 22 and local oscillator 26 is chosen; however, 
the phase relationship of the laser diodes is maintained. 

25 Different types of phased arrays are designed for 

optimal examination and imaging of different human organs 
(e.g., human head or breast). The amplitude and phase of 
the signals can be monitored on a precision oscilloscope. 
In order to scan the phased array past a fixed object of 

30 approximately known position, as in needle localization 
procedures, the location of the input and output ports 
will be scanned past the object and the position of 
maximum phase shift will be recorded in one-dimension; 
however, detection in two and three dimension can be 

35 performed in the same way. 
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In the preferred mode of operation, the array of 
sources is phased 180° apart, as shown in Fig. 8A. There 
is a sharp 180° transition of photon density wave, a 
large gradient in photon density, from S 2 , S 2 sources to 
5 the S 3 , S 4 sources. Thus, the radiated field gives an 
amplitude null and a phase transition of 180° 
corresponding to the y-z plane, i.e., perpendicular to 
the detector. If a larger number of similarly phased 
sources is used, the transitions are even sharper. The 

10 array produces a uniform photon density pattern on each 
side of the array, as shown in Figs. 8B and 8C. If an 
absorbing object is placed in this directional field of 
diffusing optical waves, imbalance in the photon density 
is measured. The detection of a hidden object is 

15 accomplished by translating the experimental transmitter- 
receiver system of Fig. 8A. 

In addition to the mechanical scanning achieved by 
moving of the input-output port system, electronic 
scanning can be performed using the multiple source and 

20 multiple detector system of Fig. 2. As shown in Fig. 2B 
for an array of five sources, there is a 180° phase 
transition in the resulting migrating field due to the 
180° phase difference between the antiphased sources 
radiating amplitude modulated light. The plane of the 

25 180° phase transition can be shifted in parallel by 
appropriately varying the apportionment of 0° and 180° 
phases on the sources. This is performed by sequentially 
switching the phase of the sources by 180°. In each 
case, the detection port located on this plane is used 

30 for collecting the data. As the sources are 

electronically switched by 180°, the detection array can 
be also electronically switched from one detection port 
to another. The signal from the receiving optical fiber 
is coupled to one shared PMT detector. However, the 

35 system can also include several detectors. If the 
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systems of Figures 1 or 1A are used, the electronic 
source scanning can be combined with synchronous 
mechanical movement of the detection port. 

In general, the invention utilizes the photon 
density gradient created in the migrating field since it 
increases the resolution of the detection. As known to 
one skilled in the art, the photon density gradient 
formed by interference effects of the individual wave can 
be created not only by appropriate phasing of the sources 
but also by other methods such as appropriately spacing 
the sources, creating an imbalance in the radiated 
intensity of the individual sources, and other. 

Fig. 8A shows the arrangement of the input ports 

11 to 17 and detection port 19 of Fig. 1. As described 

above, light of each laser diode 12 through 18 is 

intensity modulated at the 200 MHz frequency. Wavelength 

of the intensity modulated radiation is 

X = / 47rc/n f* 
\ 3f/i fl J 

wherein f is the modulation frequency of 200 MHz, /x 8 is 
the scattering factor which is approximately 10 cm" 1 in an 
Intralipid* solution with refractive index n, and c is 
3 x 10 8 cm/s. Thus, the expected wavelength is about 7 
cm. 

The input ports S 1# S 2 , S 3 , and S 4 are set 3.5 cm apart 
and are anti-phased by 180° using a push pull 
transformer. The antiphased array creates a large 
gradient in photon density chosen to take advantage of 
the destructive interference with the null detection. 
The laser diodes emitting 754 nm light are intensity 
modulated at 200 MHz using master oscillator 22, and the 
local oscillator 26 is operating at 200.025 MHz to 
perform the dynode modulation of PMT detector 24. The 
detected intensities and phase shifts of an x-direction 
scan (Fig 8A) of detection port 19 are plotted in Figs. 
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8B and 8C, respectively. As expected, the intensity has 
a sharp minimum in between sources S 2 and S 3 where the 
phase is changed 180°. The peak width at half maximum is 
about 2 cm. In addition to the x-direction scan of the 
5 detection port, the detection port was scanned in y- 
direction wherein, as expected, no variation was 
observed. 

Referring to Fig 9A, cylindrical objects of 
different diameter, d, were scanned using the previously 

10 described phased array. The objects were placed in the 
middle of the linear array displaced 2.5 cm from the x- 
axis. The detection port was located on the x-axis and 
each object was moved parallel to the x-axis at the 2.5 
cm y displacement. The intensity and phase shift 

is detected at different locations are plotted in Figs. 9B 
and 9C, respectively. The intensity pattern for each 
moving object has two maximum and one minimum when the 
scanned object was located at x=0, y=2 . 5 point during its 
scan along the x-axis. At this point, a large phase 

20 change is detected, as shown in Fig. 9C. The phase 

detection has inherently larger resolution of a localized 
absorber; a hidden object of size as small as 0.8 mm can 
be detected. 

The response due to different absorption of the 
25 hidden object was studied using a 5 mm cylinder of 

different absorption coefficient scanned by the 4 element 
phased array of Fig. 9A. The detected phase change is 
shown in Fig 9D. The 5 mm black rod displays the largest 
phase change due to its high absorption, and the cylinder 
30 filled with cardiogreen 3.5 mg/1 which has absorption 
coefficient jx a =2 00 cm" 1 shows the smallest phase change. 
In scanning of a hidden object, these experiments 
correspond to mechanically displacing the source detector 
system, or electronically scanning the subject. 
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Scanning of two objects of a different diameter is 
shown in Fig. 10A. Two cylinders of different diameter 
are scanned across the four element phased array located 
on the x-axis. The detection port in located at y = 5 
5 cm. In Fig 10B the detected phase change in plotted 
against the displacement of these objects. Curve A 
represents the phase change of two cylinders of diameters 
5mm and 10mm separated 3cm apart. Curve B was measured 
using 16mm cylinder instead the 5mm cylinder. In this 

10 case, wherein the two cylinder separation is smaller, the 
phase detector can not resolve the two objects. 

The imaging resolution is increased by increasing 
the number of elements of the phased array, since the 
main lobe of the resultant beam becomes much sharper, the 

15 gradient of photon density is larger. Phased arrays of 
different number of elements and different shapes are 
used for imaging different organs. For example, in tumor 
imaging, the four element phased array of Fig. 8A having 
an approximately linear shape can be used for imaging of 

20 the brain. On the other hand, a rectangular or a 
circular phased array would be used for imaging of a 
hidden tumor in the breast. The modulation frequency and 
the element spacing is adjusted to obtain proper 

focussing in each case. 

25 Alternative embodiments of suitable optical 

imagers are disclosed in a PCT application 
PCT/US93/05868, filed June 17, 1993, and published as 
International Publication No. WO 93/25145, which is 
incorporated by reference as if fully set forth herein. 

30 In another embodiment, the present invention 

envisions imaging systems able to calculate the average 
migration pathlengths. Referring to Figs. 1 and 1A, in 
one mode of operation, the signal from master oscillator 
22 is mixed with a set of four local oscillators 

35 operating at offset frequencies of 25, 35, 45, and 55 kHz 
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(not shown in Figs. 1 and 1A) ; there is one local 
oscillator operating at an offset frequency associated 
with each laser diode. Thus, the output of each laser 
diode is intensity modulated at the master oscillator 
5 frequency plus the frequency of its local oscillator. 
The intensity modulated radiation of each laser diode is 
simultaneously coupled to the tissue. 

Detection of the optical field is performed in the 
same way as described for the other embodiments. The 

10 detected signal is heterodyne mixed directly at the PMT 
detector. The detector outputs signals at four different 
offset frequencies associated with each diode. These 
signals are fed into the phase detector wherein the phase 
and the intensity of the detected radiation are measured. 

15 There are either four phase detectors (only one detector 
is shown in Fig. 1) operating alternatively at different 
frequencies or one phased detector is used in a time 
shared mode of operation. The phase shift and the 
intensity of a detected heterodyned signal depend on the 

20 tissue through which the scattered and absorbed radiation 
migrated. When using several radiation sources of 
selected carrier frequency and phase, the resulting 
radiation has directional properties and the detected 
intensity and phase shift depend on the pathlength along 

25 which the radiation was scattered and absorbed. The 

tissue properties are determined from the detected phase 
shift and intensity values and from the known input ports 
and detection port geometries. The measured average 
pathlengths, <L>, can also be determined. The detected 

30 phase shift is converted to an effective migration 

pathlength <L> by using the low frequency approximation 
6=27rf <L>n/c, wherein f is the modulation frequency, c is 
the speed of light (3 x 10 8 cm/s) , and n is the 
refractive index of the medium. (For detailed discussion 

35 see Analytical Biochemistry, Vol. 195, pages 330-351, 
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1991 which is incorporated by reference as if fully set 

forth herein.) 

To illustrate imaging by detecting migration 
pathlengths, we use an example of photon migration in a 
5 tissue with a strongly absorbing object, a perfect 

absorber (M a "* 00 ) of radius R. Referring to Figs. 5A, 5B, 
and 5C the distribution of pathlengths defines an optical 
field that exists between a point detector, D, and 
source, S, separated by distance p and located on the 

10 exterior of an examined tissue which is a semi-infinite, 
strongly scattering medium. As shown in Fig 5A, 
infinitely far away from the field, a perfect absorber 
does not alter the banana-shaped optical field of photons 
emitted by source S and detected at detector D. As the 

is object enters the optical field (Fig. 5B) , the photons 
which have migrated the farthest distance from D and S 
are eliminated by the absorption process inside the 
perfect absorber of radius R. Since photons which travel 
the longest pathlengths are absorbed, the approach of an 

20 object shortens the distribution of pathlengths, or 

alternatively, shortens the average pathlength <L>. As 
the object moves closer, and the optical field surrounds 
the object (Fig. 5C) , some of the detected photons have 
travelled "around" the object, which is detected as 

25 lengthening the distribution of pathlengths. Thus, the 
average pathlength measurement can reveal location of a 
strongly absorbing component of a tissue (e.g., tumor or 
localized bleeding) . 

Even though this pathlength computation approach 

30 requires in most cases extensive computational 

capabilities, it can yield useful information in the 
localization procedures and can provide an useful 
supplement to the above described directional approach. 

In another preferred embodiment, shown in Fig. 11, 

35 a low frequency imaging system 100 includes two light 
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sources 102 and 104 (e.g., tungsten lamps, LEDs) of 
electromagnetic radiation of a visible or infrared 
wavelength. The sources 102 and 104 introduce into 
examined tissue 8 at input locations 101 and 103, 
5 respectively, the radiation of known intensities that 
define a null plane 105 in tissue 8. A detection port, 
optically connected to a detector 106, is located at a 
detection location 107 of null plane 105. Detector 106 
detects sequentially radiation that has migrated from the 

10 input locations 101 and 103 to detection location 107. A 
detector circuit 110, connected to detector 106, receives 
sequentially two detection signals 108, the first 
corresponding to radiation that has migrated from first 
input location 101 to detection location 107 and the 

15 second corresponding to radiation that has migrated from 
second location 103 to detection location 107. Detector 
circuit 110, includes a sample-and-hold circuit and a 
subtraction circuit constructed to subtract the first 
detection signal from the second detection signal. A 

20 processor 120 controls the entire operation, receives the 
differential signal 119 from detector circuit 110, 
governs the operation of intensity controller 13 0, and 
sends imaging data to a display 140. 

Intensity controller 130 regulates the timing and 

25 the intensity on the radiation emitted from each source. 
For equal intensities of the emitted radiation, null 
plane 105 is located at a midpoint between the two 
optical field patterns generated from sources 102 and 
104 . Detector 106 is located at the null plane or at a 

30 known location relative to the null plane. 

Referring to Fig. 11A, detector circuit 110 
receives detection signal 108 from a diode detector 106. 
The circuit enables correction for the dark current/ noise 
that comprises background light, DC offset of the 

35 operational amplifiers, photodiode dark current, 
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temperature effects on the outputs of individual 
components and variations due to changing environment. 
The circuit # also enables subtraction of the first 
detection signal from the second detection signal. 
5 The system performs data acquisition in two 

cycles , each having four steps synchronized by its 
internal oscillator. In the first cycle, detecting light 
that has migrated from input location 101 to detection 
location 107, the first step is performed by having the 

10 light sources off. The output is directed to an 

integrator 112c and integration capacitor 113c is charged 
to the dark level voltage. In the second step, the first 
light source is turned on- The preamplifier output that 
corresponds to the intensity of the detected light is 

15 directed to integrator 112c in a way to charge capacitor 
113c with current of polarity opposite to the polarity of 
the charging current in the first step. This is achieved 
using appropriate ON/OFF combination of switches A and B. 
The voltage of capacitor 113c is charging to a value 

20 which, at the end of this step, represents the total 

signal minus the dark level noise signal. In the third 
step, both switches A and B are turned OFF to disconnect 
both the positive unity gain and the negative unity gain 
operational amplifiers (112a and 112b). Then, the output 

25 of integrator 112c is moved via switch C to a hold 

circuit with integrator 112d which also functions as a 
low pass filter. The output is the first detection 
signal corrected for the background noise. In the fourth 
step, the switches A, B and C are open and switch D is 

30 closed in order to discharge capacitor 113c through a 47K 
resistor. At this point, the circuit of integrator 112c 
is reset to zero and ready for the first step of the 
second detection cycle. 

The second detection cycle is again performed in 

35 four steps employing the second light source (L2) instead 
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of the first light source (LI) . After the second 
detection signal, corrected for the background noise, is 
obtained, the two detection signals are subtracted in 
module HOB and the deferential signal 119 is sent to 
5 processor 120. 

Spectrophotometer 100 has two equivalent 
arrangements. The first arrangement shown in Fig. 11 
includes two sources 102, 104 (or a single source 
optically coupled to two input ports) and single detector 

10 106 with detection port 107 located on the null plane 

105. The second arrangement includes a single source and 
two detectors. The source is optically coupled to an 
input port located on the null plane defined by the 
positions of the two detectors (or a single detector 

15 optically coupled to two detection ports) . For 

substantially "homogeneous" tissue (i.e., normal tissue) 
or if a tumor is located symmetrically on the null plane, 
the detectors will detect a symmetric signal. Otherwise, 
the detectors will detect a asymmetric signal due to the 

20 modification of the optical field caused by the tumor 

having different scattering or absorption properties than 
the normal tissue. 

The source-detector geometry may be a transmission 
geometry as shown in Fig 11 or a reflection geometry as 

25 shown in Fig. 14A. Another source-detector probe that 
utilizes a reflection geometry is disclosed in 
International Publication No. WO 92/20273, filed May 18, 
1992, which is incorporated by reference as if fully set 
forth herein. 

30 To examine a volume of tissue, the system may 

sweeps the null plane position over the examined tissue 
by changing the relative value of the intensities of the 
emitted radiation. Detector 106, which may include an 
array of detectors, is again located at the null plane or 

35 at a known location relative to the null plane. 
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Alternatively, the null plane scanning is achieved by 
moving the input ports and the detection port to 
different locations, or by scanning the introduced light 
over the tissue surface by orienting a mirror. 
5 Figs. 12 and 13 display another embodiment of the 

low frequency imaging system that utilizes a two 
dimensional array of light sources. Referring to Fig. 
12, imaging system 140 includes four sources 142, 144, 
146, and 148 (labeled W, N, E, and S) which are turned on 

10 and off as shown in a timing diagram 150. The system 
uses two differential detectors 152 and 154 , located 
relative to the null plane, which produce deferential 
signals 119 and 199a, respectively, corresponding to the 
radiation emitted from sources 144, 148, and sources 142, 

15 146. Referring to Fig. 13, imaging system 150 utilizes a 
two dimensional, dual wavelength array. The system 
operates based on similar principles as described for the 
imaging systems of Figs. 11 and 12. 

To increase resolution of the system an exogenous 

20 pigment sensitive to the introduced radiation (a light 
sensitive contrast agent) can be introduced, e.g., by 
injection, into the examined tissue. Then, the imaging 
system uses radiation of a wavelength sensitive to the 
exogenous pigment that also may be preferentially 

25 absorbed in a tumor located in the examined tissue. 
Alternatively, the imaging system uses radiation of a 
wavelength sensitive to an endogenous pigment of the 
examined tissue. 

Figs. 14 and 14A show a scanning system 160 

30 constructed for imaging of breast tissue. Scanning 

system 160 may employ spectroscopic systems of Figs. 1, 2 
or 3, or systems of Figs. 11, 12 or 13. Scanning system 
160 includes an optical coupler 162, which may have 
cubical or cylindrical shape and is filled with optical 

35 medium 164. Optical coupler 162 is positioned over the 
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breast: near the chest wall. As described above in U.S. 
Patent 5,402,778, which is incorporated by reference, the 
optical properties, pressure and volume of medium 164 may 
be controlled by an external system connected to the 
5 coupler by a set of tubes. The optical matching fluid 
(e.g., twice-diluted J&J baby lotion) is contained within 
pliable, optically an transparent barrier. The inside 
walls of coupler 162 may be coated with a film that 
reflects light in the visible or near infra-red range 

10 back to the matching fluid to prevent escape of photons 
from the tissue surface. The optical coupler may have 
different sizes or may have an adjustable volume so that 
the coupler can have a selected distance between the 
breast surface and the inside walls. (The preferred 

15 distance is about 1 centimeter, but for a very small 
tissue a larger distance is preferable to achieve semi- 
infinite boundary conditions.) Thus the coupler is also 
useful for examination of the breast of a small size or 
after a surgical removal of the breast tissue. After 

20 placement of coupler 162, the volume of medium 164 is 
adjusted so that the barrier fits snugly around the 
examined breast. Alternatively, the optical medium is a 
pliable solid, for example, an absorbing gel containing 
metallic or oxide spherical particles, silky glass beads 

25 as scatterers or a suitable plastic material. 

Fig. 14 A depicts a set of couplers 162A and 162B 
for simultaneous scanning of both breasts. Attached to 
each coupler are source-detector probes (168A, 168B, 
168C, 168D, 169A, 169B, 169C, 168D) , which include one or 

30 more optical sources or detectors described above. The 
probes are movable on a rail 170. In an automatic 
positioning system, each probe is connected to a servo 
motor (step motor) that is operated by a controller. 
Depending on the spectroscopic system, a fiber 172 may be 

35 used to collect, at a detection port 174, radiation that 
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has migrated in the examined tissue and couple the 
radiation to a detector. Alternatively, fiber 172 may be 
used to couple, at input port 174, radiation to the 

examined tissue. 
5 In an electro-optic scan, a computer controller 

maintains coordinated positions of the probes to the 
selected combination of the transmitters and receivers. 
The scan is performed on a single breast or 
simultaneously on the contralateral breast. The 

10 sensitivity of the simultaneous scan is increased be 

measuring a differential signal. A computer displays the 
detected signal or the differential signal in a 3 
dimensional coordinate system. To increase the 
resolution, a contrast agent (e.g., cardio-green, 

15 indocyanine-green) which is preferentially accumulated in 
a tumor may by injected intravenously. Several scans are 
performed to observe the time dependence of the decay and 
identify a location of a suspected anomaly. The system 
can also calculate the scattering coefficient and 

20 absorption coefficient of the suspected anomaly as in 
U.S. Patent 5,402,778 and the reference cited therein. 

The scan may be combined with a needle 
localization procedure, X-ray mammography or an MRI scan 
as is described in International Application WO 95/02987 

25 (PCT/US94/07984, filed July 15, 1994). 

What is claimed is: 
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1. A spectroscopic method of deriving data 
signals useful for examination of biological tissue of a 
subject, comprising the steps of: 

providing at least one source of electromagnetic 
5 radiation of a visible or infrared wavelength selected to 
be scattered and absorbed while migrating in the tissue, 
said source being optically coupled to at least two input 
ports, and providing a detection port optically coupled 
to a detector, said detector connected to a detector 
10 circuit; 

positioning a first input port and a second input 
port relative to selected input locations of a subject; 

selecting for each input port first and second 
radiation intensities to be introduced to the tissue, 
15 said selected radiation intensities defining a null plane 
in the tissue; 

positioning said detection port relative to a 
selected detection location of the examined tissue 
corresponding to said null plane, said input locations 
20 and detection location defining a volume of the examined 
tissue of the subject; 

introducing into the tissue, at said first input 
port, radiation of said first intensity; 

detecting, at said detection port, said first 
25 radiation that has migrated in the examined tissue; 

storing, in said detector circuit, a first 
detector signal corresponding to said first detected 
radiation; 

introducing into the tissue, at said second input 
30 port, radiation of said second intensity; 

detecting, at said detection port, said second 
radiation that has migrated in the examined tissue; 

storing, in said detector circuit, a second 
detector signal corresponding to said second detected 
35 radiation; and 
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subtracting said first detector signal from said 
second detector signal to obtain processed data 
corresponding to properties of said tissue volume. 

2. A spectroscopic method of deriving data 
5 signals useful for examination of biological tissue of a 
subject, comprising the steps of: 

providing a source of electromagnetic radiation of 
a visible or infrared wavelength selected to be scattered 
and absorbed while migrating in the tissue, said source 
10 being optically coupled to an input port, and providing 
at least two detection ports optically coupled to at 
least one detector, said detector connected to a detector 
circuit; 

positioning said input port relative to selected 
15 input locations of the tissue; 

positioning a first detection port and a first 
detection port relative to selected detection locations 
of the examined tissue, said locations defining a null 
plane in the tissue a volume of the examined tissue of 

20 the subject; 

introducing into the tissue, at said input port, 

radiation of a selected intensity; 

detecting, at said first detection port, radiation 
that has migrated in the examined tissue and storing, in 
25 said detector circuit, a first detector signal 
corresponding to said detected radiation; 

detecting, at said second detection port, 
radiation that has migrated in the examined tissue and 
storing, in said detector circuit, a second detector 
30 signal corresponding to said detected radiation; and 

subtracting said first detector signal from said 
second detector signal to obtain processed data 
corresponding to properties of said tissue volume. 
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3 • The spectroscopic method of claim 1 or 2 
wherein said first detector signal and said second 
detector signal are stored and subtracted in an analog 
form by said detection circuit. 

5 4. The spectroscopic method of claim 1 or 2 

further comprising, before said storing steps, converting 
said first and second detection signals to a digital 
form, said subtracting step being performed digitally on 
said first detector signal and said second detector 
10 signal in an digital circuit. 

5. The spectroscopic method according to any 
preceding claim wherein said steps of positioning said 
input and detection port at the respective selected 
locations include placing said ports on the surface of 

15 the examined tissue. 

6. The spectroscopic method of claim 1 or 2 
wherein said steps of positioning said input and 
detection ports include orienting said ports relative to 
said input and detection locations. 

20 7. The spectroscopic method of claim 5 wherein 

said orienting said input ports include directing an 
optical element of said port to said location. 

8. The spectroscopic method according to any 
preceding claim wherein said positioning steps are 

25 performed in a manner that sweeps said null plane and 
said detector over at least a portion of said volume of 
the examined tissue. 

9. The spectroscopic method of claim 1 wherein 
said step of selecting said first and second radiation 
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intensities is performed in a manner that sweeps said 
null plane over at least a portion of said volume of the 
examined tissue and said method further comprising , 
simultaneously with said sweeping, positioning said 
5 detection port to detection locations corresponding to 
said swept null plane. 

10. The spectroscopic method of claim 2 wherein 
said detecting steps are performed in a manner that 
changes the relative sensitivity of detection at said 

10 first and second detection port in order to sweep said 
null plane over at least a portion of said volume of the 
examined tissue and said method further comprising, 
simultaneously with said sweeping, positioning said input 
port to input locations corresponding to said swept null 

15 plane. 

11. The spectroscopic method according to any 
preceding claim further comprising locating, in said 
tissue volume, a tissue region exhibiting different 
scattering or absorptive properties than the rest of the 

20 examined tissue volume. 

12. The spectroscopic method of claim 11 further 
comprising imaging the examined tissue including said 
tissue region of different scattering or absorptive 
properties. 

25 13 . The spectroscopic method of claim 1 or 2 

further comprising displaying an image of the examined 
tissue by utilizing to said processed data and relative 
locations of said input ports and said output ports. 

14 . The spectroscopic method according to any 
30 preceding claim further comprising a step of introducing 
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an exogenous pigment: into the tissue, said selected 
wavelength being sensitive to said pigment. 

15. The spectroscopic method of claim 14 wherein 
said exogenous pigment is preferentially accumulated in a 
tissue region exhibiting different scattering or 
absorptive properties. 

16. The spectroscopic method of claim 14 wherein 
said exogenous pigment is a fluorescing pigment when 
irradiated by a selected wavelength, and said detecting 
steps detecting preferentially radiation wavelength of 
said fluorescing pigment. 
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17 . A spectroscopic system for examination of 
biological tissue of a subject, comprising: 

at least one light source of electromagnetic 
radiation of a visible or infrared wavelength; 
5 at least two input ports, optically coupled to 

said light source, constructed to introduce at selected 
input locations of the examined tissue said radiation of 
known intensities that define a null plane in the tissue; 

a detection port located at a selected detection 
10 location of the examined tissue relative to said null 
plane; 

a detector, optically coupled to said detection 
port, constructed to detect radiation that has migrated 

in the examined tissue; 
is a detector circuit connected to and receiving 

detection signal from said detector, said detector 

circuit including a sample-and-hold circuit and a 

subtraction circuit; 

said detector circuit constructed to correlate 
20 emission of said radiation from a first input port with 

detection of radiation scattered and absorbed while 

migrating in the tissue, said detected radiation being 

stored as a first detection signal; 

said detector circuit further constructed to 
25 correlate emission of said radiation from a second input 

port with detection of radiation scattered and absorbed 

while migrating in the tissue, said detected radiation 

being stored as a second detection signal; 

said subtraction circuit constructed to subtract 
30 said detection signals; and 

a processor, connected to and receiving said 

processed data from said subtraction circuit, constructed 

to evaluate the examined tissue. 
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18. A spectroscopic system for examination of 
biological tissue of a subject, comprising: 

a source of electromagnetic radiation of a visible 
or infrared wavelength; 
5 an input port, optically coupled to said light 

source, constructed to introduce at a selected input 
location of the examined tissue said radiation; 

a detector optically coupled to at least two 
detection ports located at selected detection locations 
10 defining a null plane in the examined tissue, said 
detector constructed to detect radiation that has 
migrated in the examined tissue to said detection ports; 

a detector circuit connected to and receiving 
detection signal from said detector, said detector 
15 circuit including a sample-and-hold circuit and a 
subtraction circuit; 

said detector circuit constructed to correlate 
emission of said radiation from said input port with 
detection of radiation scattered and absorbed while 
20 migrating in the tissue at said first detection port, 

said detected radiation being stored as a first detection 
signal; said detector circuit further constructed to 

correlate emission of said radiation from said input port 
with detection of radiation scattered and absorbed while 
25 migrating in the tissue at said second detection port, 
said detected radiation being stored as a second 
detection signal; 

said subtraction circuit constructed to subtract 

said detection signals; and 
30 a processor, connected to and receiving said 

processed data from said subtraction circuit, constructed 
to evaluate the examined tissue. 



19. A spectroscopic system of claim 17 further 
comprising intensity control means constructed to 
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regulate intensities of radiation emitted from said first 
and second input ports. 

20. A spectroscopic system of claim 19 wherein 
said intensity control means is further constructed to 
5 regulate said intensities in a manner that sweeps said 
null plane over at least a portion of said volume of the 



21. A spectroscopic system of claim 19 or 20 
further comprising positioning means constructed to 

10 displace said detection port to detection locations 
corresponding to said null plane. 

22. A spectroscopic system of claim 17 , 19, 20 or 
21 further comprising positioning means constructed to 
displace said input ports to selected locations. 

15 23. A spectroscopic system of claim 18 further 

comprising detector controller means constructed to 
changes the relative sensitivity of detection at said 
first and second detection port in order to sweep said 
null plane over at least a portion of said volume of the 

20 examined tissue. 

24. A spectroscopic system of claim 18 or 23 
further comprising positioning means constructed to 
displace said input port to a selected location. 



25. A spectroscopic system according to any 
25 preceding system claim wherein said input or detection 
ports are arranged in a linear array. 
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26. A spectroscopic system according to any 
preceding system claim wherein said input or detection 
ports are arranged a two dimensional array. 

27. A spectroscopic system according to any 

5 preceding system claim wherein said subtraction circuit 
further includes an analog to digital converter, 
connected to said sample-and-hold circuit, constructed to 
digitize said first and second detection signals to 
produce digital detection signals, said subtraction 
10 circuit subtracting said digital detection signals to 
obtain said processed data. 

28. A spectroscopic system according to any 
preceding system claim further comprising 

an image processor connected to receive said 
15 processed data from said processor and create image data; 
and 

a display, connected to said image processor, 
constructed to display said image data representing the 



20 29. A spectroscopic system according to any 

preceding system claim wherein said processor further 
constructed to locate, in said tissue volume, a tissue 
region exhibiting different scattering or absorptive 
properties than the rest of the examined tissue volume. 

25 3 0. A spectroscopic system of claim 29 wherein 

said image processor is further constructed to create 
image data of said located tissue region exhibiting 
different scattering or absorptive properties; and 

said display constructed to display said image 

30 data . 
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31. A spectroscopic system according to any 
preceding system claim wherein said wavelength is 
sensitive to an endogenous pigment of the examined 
tissue. 

5 32. A spectroscopic system according to any 

preceding system claim wherein said wavelength is 
sensitive to an exogenous pigment of the examined tissue. 
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